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ABSTRACT 

Based on an ultra deep (230 ks) ROSAT HRI imaging of M101, we have 
detected 5 X-ray sources that coincide spatially with optical emission line 
features previously classified as supernova remnants in this nearby galaxy. 
Two of these coincidences (SNR MF83 and NGC5471B) most likely represent 
the true physical association of X-ray emission with shock-heated interstellar 
gas. MF83, with a radius of ~ 134 pc, is one of the largest remnants known. 
NGC5471B, with a radius of 30 pc and a velocity of at least 350 km s _1 (FWZI), 
is extremely bright in both radio and optical. The X-ray luminosities of these 
two shell-like remnants are ~ 1 and 3 x 10 38 ergs s _1 (0.5-2 keV), about an 
order of magnitude brighter than the brightest supernova remnants known 
in our Galaxy and in the Magellanic Clouds. The inferred blastwave energy 
is ~ 3 x 10 52 ergs for NGC5471B and ~ 3 x 10 53 ergs for MF83. Therefore, 
the remnants likely originate in hypernovae, which are a factor of ^ 10 more 
energetic than canonical supernovae and are postulated as being responsible for 
7-ray bursts observed at cosmological distances. The study of such hypernova 
remnants in nearby galaxies has the potential to provide important constraints 
on the progenitor type, rate, energetics, and beaming effect of 7-ray bursts. 



Subject headings: galaxies: general — galaxies: individual (M101) — ISM: 
supernova remnants — X-rays: general — gamma-ray: bursts 



1. Introduction 

The face-on grand-design spiral M101 (Fig. 1) is an ideal site for studying various 
X-ray source populations and their relationship to other galactic components. At a 
well-determined distance D = 7.2 ± 0.4 Mpc (thus T = 2.1 kpc; Stetson et al. 1998), 
the galaxy is probably the largest in the sky, except for members of our Local Group. 
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The galaxy is also in a direction of exceptionally low foreground Galactic extinction 
[E(B — V) ~ 0; N H ~ 1.1 x 10 20 cm -2 ]. A recent optical survey of the galaxy by Matonick 
& Fesen (1997, MF hereafter) has identified 93 supernova remnant (SNR) candidates. The 
galaxy also contains numerous giant HII regions (GHRs) and complexes (GHCs), located 
primarily in outer spiral arms. Some of these GHRs are substantially brighter than the 
LMC 30 Dor nebula, which is the brightest HII region in the Local Group. Thus, M101 
provides a particularly favorable circumstance to explore starburst-related phenomena. 

Based on an ultra-deep image from the ROSAT High Resolution Imager (RHRI), 
Wang et al. (1999) have detected 51 sources with X-ray luminosities in the range of 0.4 
- 20 xlO 38 ergs s _1 (0.5-2 keV). A cross-correlation of these X-ray sources with SNR 
candidates in the galaxy leads to a surprising discovery of 5 position coincidences. Unlike 
several previously-known X-ray-bright SNRs (in other nearby galaxies), which result from 
strong interaction between SN ejecta with dense circumstellar media, these apparently 
X-ray-emitting remnants in M101 are optically-resolved, and are thus relatively well evolved. 
Most importantly, the observed X-ray emission, if arising in shock-heated interstellar gas, 
implies the blastwave energy of each remnant is a factor of ^ 10 greater than the energy 
(~ 10 51 ergs; e.g., Hughes et al. 1998) expected for an SN. Therefore, the remnants may 
well be the products of hypernovae, as have been proposed to explain 7-ray bursts (GRBs) 
and their afterglows (Paczyhski 1998; Fryer & Woosley 1998). In this Letter, we describe 
this discovery and explore implications of such remnants in the study of GRBs. The term 
"hypernova", proposed first by Paczyhski (1998) for a super-energetic GRB/afterglow event, 
is loosely used here to refer an explosion that appears too energetic to be a nominal SN. 



2. Cross-Correlation between X-ray Sources and SNR Candidates 

Table 1 lists 5 position coincidences between X-ray sources and known SNR candidates 
in M101 (see also Fig. 1). Except for NGC5461B, which was discovered in radio continuum 
(Skillman 1985; Yang, Skillman, & Sramek 1994) and was confirmed in optical (Chu & 
Kennicutt 1986; Chen & Chu 1999), all other SNRs used in the cross-correlation are from 
the MF survey, which consists of both optical imaging and spectroscopic observations. The 
imaging observations covered five partially-overlapped regions of 8f7x8'7 each. Emission 
nebulae with [S II]/Ho; ^ 0.45 were identified as SNRs. Spectroscopic followup observations 
were carried out for 14 objects, including MF83 in Table 1. The SNR positions (J2000; 
Table 1) are accurate to 1" , radii R to about ±5 pc, and both Ha Luminosities L(H a ) 
and [S II] /Ha ratios to 15% (MF). The 51 X-ray sources were detected with signal-to-noise 
ratios > 3.5 within a field of 17' radius covered by the RHRI image (Wang et al. 1999). 
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The RHRI was sensitive to photons in the 0.1-2 keV range and had a point spread function 
(PSF) with a half power-encircled radius ranging from ~ 2'/5 on-axis to 20" at 17' off-axis. 
We corrected the astrometry of the image by matching centroids of 7 point-like X-ray 
sources with their optical counterparts. Table 1 includes only the SNRs whose position 
offsets A-r-o from X-ray source centroids are less than twice the la uncertainty radii r err 
(Wang et al. 1999). Table 1 also includes the projected galactocentric radius R g , expanding 
velocity v e , and optical morphology as well as the X-ray source number and 0.5-2 keV 
luminosity of each remnant. 

What is the probability for a chance superposition? The MF survey had a complete 
coverage of the central field of R g ^ 8' from the M101 nucleus, which is at R.A., Dec. 
(J2000) = 54 A 3 m 12 i !7, 54°20'54" (Israel, et al. 1975). This field includes 29 X-ray sources 
and 87 SNRs. Our adopted 2a matching radii enclose a total region of 583 arcsec 2 . The 
statistical probability to find by chance 1 match in the region, for example, is 0.45. Four 
SNRs in Table 1 are within the field. It is quite possible that one or two of the coincidences 
in Table 1 are just the chance superposition. Indeed, H19 shows strong time variability 
and has a relatively hard spectrum (Wang et al. 1999). Thus the source most likely 
represents an accreting X-ray binary system. Although this system may still be related to 
the explosion responsible for the remnant, the X-ray emission is not due to shock- heated 
interstellar gas. To find 4 random SNR/X-ray source superpositions in the field, however, 
the probability falls to ~ 10~ 3 . Thus, we expect that at least some of the coincidences in 
Table 1 do represent physical associations between the X-ray sources and SNRs. 

The most likely true X-ray-bright remnants in Table 1 are MF83 and NGC5471B. 
Their A x _ are the smallest, well within the la error radii. Both of them are shell-like (the 
other three have filled morphology; Table 1). Most importantly, these two SNRs are known 
to be the most unusual in M101 (MF; Yang, Skillman, & Sramek 1994; Chen & Chu 1999). 
MF83 is one of the largest SNRs known. Its high [S II] /Ha ratio and bright [OIII] emission 
clearly suggest that the optical nebula is shock heated; the [O III]/H J g(= 2.0) ratio indicates 
a shock velocity of ~ 86 km s _1 (Table 17 in MF). The remnant, located in a region just 
off a spiral arm, shows no evidence for an association with any concentration of OB stars. 
NGC5471B is an extended, nonthermal radio source with a flux of 4.4±0.8 mJy at 30 cm 
and a spectral slope of -0.73 (Yang, Skillman, & Sramek 1994). Based on optical echelle 
spectroscopy, Chu & Kennicutt (1986) detected shocked gas with velocity ^ 350 km s _1 
(FWZI) and inferred the required supernova blastwave energy as ~ 10 52 ergs. Therefore, 
the identification of both remnants are firm and they all appear to be unusually energetic. 
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3. Constraints on Remnant Parameters 

To characterize the spectral properties of X-ray-emitting plasma, we have further 
analyzed a ROSAT PSPC observation of M101 (Wang et al. 1999). Though only with 
a limited spatial resolution of ~ 30" FWHM on-axis, this observation provides a spectral 
resolution of 5E/E ~ 0.43(0.93 keV/£)°- 5 . The PSPC spectral data for NGC5471B is 
reasonably reliable, as it is located in a relatively isolated region, 22 kpc off the galaxy's 
nucleus. We find that the Raymond & Smith thermal plasma model provides a satisfactory 
fit to the data (x 2 /d.o.f. = 6.2/11), giving the plasma temperature as 0.29(0.22-0.45) and 
the absorbing-gas column density as 1.8(0.3 — 6.1) x 10 20 cm -2 ; the parameter intervals are 
all at the 90% confidence. The solar metal abundance is assumed for X-ray-absorbing gas, 
which affects only the absolute value of the column density. The normalization of the model 
is 

K = [10- 14 /(4tt j D 2 )] J n 2 e dV ~ 8.2 x 10~ 5 cm" 5 , (1) 

where n e is the electron density of the plasma, and the integration is over the entire 
volume of the remnant. The metal abundance ( of the X-ray-emitting plasma is assumed 
as 13% solar, as inferred from accurate measurements of oxygen abundance in NGC5447 
(Torres-Peimbert, Peimbert, & Fierro, 1989). The temperature of X-ray-emitting plasma 
indicates an expanding shock velocity of ~ 5.0 x 10 2 km s -1 , consistent with the lower limit 
inferred from the optical kinematics. This, together with the relatively small optical size of 
NGC5471B, suggests the remnant is in the Sedov- Taylor phase. 

We can thus further estimate physical parameters of NGC5471B, using the Sedov 
solution: 

J R = (20pc)(2E 52 t 2 /n ) 1/5 , (2) 

where R, E 52 , t 4 , and n are the radius (in units of pc), blastwave energy (10 52 ergs), age 
(10 4 yrs), and preshock density (cm -3 ) of the remnant. The shocked gas temperature can 
be expressed as 

kT ~ (1.8 x 10~ 3 keV)(R/U) 2 , (3) 
and the postshock-to-preshock density ratio is n e /n = 4. From (1), (2), and (3), we obtain 

n ~ 2 cnr 3 (i?/30 pc)' 3/2 (K/8 x 10~ 5 cm" 5 ) 1 / 2 , (4) 

E ~ (3 x 10 52 ergs)(i?/30 pc) 3/2 (K/8 x 10~ 5 cnT 5 ) 1/2 (A;T/0.3 keV), (5) 

and 

t~(2x 10 4 yrs)(A;T/0.3 keV)- 1/2 (i?/30 pc). (6) 

This age is smaller than the time for the pressure-driven snowplow (PDS) phase to start 
(Cioffi, McKee, & Bertschinger 1988), t PDS ~ (5 x 10 4 yrs)no^ 7 -E52 14 Cai3 14 ? consistent 
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with our assumption of the Sedov- Taylor phase. We further estimate the total thermal 
energy as ~ 10 52 ergs and the gas cooling rate as ~ 10 39 ergs s -1 . The inferred cooling 
time is ~ 3 x 10 5 yrs, and thus the remnant is apparently not in a radiative phase. The 
estimate of E, in particular, is not sensitive to the uncertainty in the spectral parameters, 
because K anti-correlates with kT (Fig. 2). Thus, NGC5471B most likely originates in a 
super-energetic explosion — a hypernova. 

Similarly, we can combine the optical and X-ray results to constrain the physical 
properties of MF83. Its large radius and relatively low shock velocity suggest 
that this remnant has evolved into the PDS phase, which starts at the radius 
Rpds = (32 pc)E^ 7 n^ 3/7 Ql /7 and velocity v PDS = (376 km s -1 )£^ 14 nJ /7 CS.3 14 . C = 0.3 
is assumed to be the same as the abundance measured for the nearby GHC NGC5461 
(Torres- Peimbert, Peimbert, & Fierro, 1989), and the parameter values to be inferred are 
not sensitive to this assumption. The PDS evolution (Cioffi, McKee, & Bertschinger 1988) 
follows 

R = RpDs[WW/tpDs) - l/3] 3/1 ° (7) 

and 

v e = v PD8 [(4/3)(t/tpDs) - l/3]- 7/1 °. (8) 

From these two equations and the values of R and v e (Table 1), we get E 52 ~ 13rtQ 2 . If 
n ^ 0.1 cm" 2 , MF83 is then a hypernova remnant (HNR) as well. We may estimate n 
from the X-ray data on MF83. However, the remnant is not bright enough in X-ray and 
is in a relatively crowded region; no reliable X-ray spectral constraint is yet available. 
Nevertheless, if the X-ray spectrum (e.g., the plasma temperature) is not too much different 
(a factor of <^ 3) from that of NGC5471B, the X-ray luminosity of MF83 suggests a total 
thermal energy E p ~ 6 x 10 52 ergs. Taking into account the adiabatic energy loss during the 
PDS phase, E ~ E T (R/R PDS ) 2 , we obtain n = (32/ R) 46 / 15 (37Q/v e ) u / 5 E^° ~ 2 cm" 3 ; 
thus, t ~ 1 x 10 6 yrs and E~ 3 x 10 53 ergs. 

While the above calculations are based on a simple SNR evolution model, the 
consideration of potential complications may significantly alter the estimates of the 
parameters. For example, if an explosion occurs in a pre-existing low-density (stellar wind) 
bubble, the actual blastwave energy can then be a factor of a few lower than our estimate 
(e.g., Hughes et al. 1998). The expanding velocity of MF86 is also uncertain. But we do 
not expect that such complications would affect our results qualitatively. The strong X-ray 
emission as well as the unusual optical and radio characteristics of NGC5471B and MF83 
all suggest that they originate in super-energetic explosions. 
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4. Hypernova Remnants and GRBs 

The detection of HNRs in nearby galaxies may have important implications for our 
understanding of GRBs. The recent discovery of afterglows strongly suggests that GRBs 
represent the most extreme high-energy events we have witnessed in the universe. In the 
extreme case (GRB 990123), the inferred isotropic energy release (^ 3 x 10 54 ergs) exceeds 
the rest mass of a neutron star (Kulkarni et al. 1999), although the beaming of the emission 
may substantially reduce the energy requirement (e.g., Meszaros, Rees, & Wijers 1998 
and references therein). Furthermore, such an event may contain 10 times more energy 
in outflowing non-relativistic matter (Loeb & Perna 1998), which is not observable during 
a GRB and its afterglow. Most plausible models of GRBs involve either the coalescence 
of two compact objects — two neutron stars or a neutron star and black hole — or 
collapses of massive stars and/or their mergers with compact companions (Paczyhski 1998; 
Fryer & Woosley 1998). The later scenarios have become particularly attractive, because 
of the evidence that GRB/afterglow events appear close to star- forming regions. In all 
these models, an event leads to the formation of a black hole and provides an extractable 
energy of ~ 10 54 ergs (Meszaros, Rees, & Wijers 1998). This naturally raises the hope 
to find in nearby galaxies relics of similar events as they must have interacted with the 
interstellar medium (ISM). Indeed, it has been suggested recently by Efremov, Elmegreen, 
& Hodge (1998) and by Loeb & Perna (1998) that some of H I supershells or holes in 
the ISM may represent such relics. As demonstrated in these papers, however, it is very 
difficult to distinguish the GRB scenario from other more conventional interpretations of 
the supershells: the combined action of stellar winds and SN of OB associations or impacts 
of high-velocity clouds. OB associations, for example, can be hard to detect within H I 
supershells with typical ages ^ 10 7 yrs. 

In contrast, a relatively young HNR, or a young GRB relic, in a nearby galaxy can 
be identified relatively easily and can be studied in great details at mult i- wavelengths. 
Except for its substantially greater blastwave energy, an HNR resembles an SNR and 
should have such distinct signatures as nonthermal radio continuum, high [S II] /Ha ratio, 
large shock velocity, and strong soft X-ray emission, as observed from NGC5471B and 
MF83. These signatures are not expected for a supershell powered even by a young OB 
association. Energetically, one can argue that NGC5471B or MF83 might be a result of 
multiple SN explosions. But this would require an explosion rate of ^ 1 per 10 4 yrs, or a 
mean mechanical energy input of ^ 10 40 ergs s _1 , within the boundaries of each remnant. 
Such concerted energy input could only occur in very massive OB associations, more 
massive than 30 Dor in the Large Magellanic Cloud (e.g., Wang 1999). To be consistent 
with the sizes of NGC5471B and MF83, the massive OB associations would also have to 
be very young (less than several times 10 6 yrs). Such young OB associations, however, 
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would produce luminous HII nebulae, similar to the 30 Dor nebula with predominantly 
thermal characteristics, which do not fit the observations of either NGC5471B or MF83. 
One could further recognize the associations in optical images. But, there is no evidence for 
their presence within both NGC5471B and MF83 (Chu 1999, privately communications). 
Similarly, one can distinguish young HNRs from small ring-like or half-shell features which 
might be produced by recent impacts of high-velocity clouds. At an expected velocity 
comparable to the galactic disk rotation, the impact could hardly heat gas to a temperature 
of a few times 10 6 K, which should be typical for young HNRs. In conclusion, one should 
be able to identify relics of GRBs if they are present as young HNRs in nearby galaxies. 

Could NGC5471B and MF83 be the relics of GRBs statistically? From our sample of 
2 HNRs with their ages within ~ 1 x 10 6 yrs, we may infer a rate of order 2 hypernovae 
per 10 6 yrs in M101. Within the obvious uncertainties, this rate agrees with theoretical 
estimates from stellar evolution models for binary neutron star mergers or induced neutron 
star collapses in binary systems for an L* galaxy (e.g. Qin et al. 1998), and may also 
be comparable to estimated hypernova rates, which depend largely on the history of 
star formation (e.g., Fryer & Woosley 1998). Observationally, GRBs appear to be about 
10 4 — 10 6 rarer than SN (Wijers et al. 1998). Assuming an SN rate of 1 per 50 yrs (MF) 
in M101, one then expects 0.02-2 GRBs per 10 6 yrs. Our inferred rate is at the upper end 
of this range. This may be explained by two important factors. First, the SN rate in M101 
is of great uncertainty; the MF survey covers only part of the galaxy and a substantial 
number of SNRs as expected in GHCs are not detected (MF). The true SN rate may be 
considerably higher in M101. Second, only a fraction of hypernovae are likely observable 
as GRBs at cosmological distances, presumably due to a moderate beaming effect of 7-ray 
emission as expected in GRB models. Therefore, the detection of the 2 remnants in M101 
is reasonably consistent with the GRB statistics. 

Clearly, our detection of the 2 apparently young HNRs in M101 represents only the 
first step to firmly establish their identity. No matter what this identity might be, however, 
these remarkable remnants must originate in very energetic events. Other nearby galaxies 
may also contain similar remnants. In particular, there is a population of persistent X-ray 
sources with luminosities greater than the Eddington limit for a ~ 1.6M neutron star. 
Such superluminous sources tend to be associated with galaxies of active massive star 
formation, and are absent in members of the Local Group (e.g., Wang et al. 1999). It is 
conceivable that some of these sources, as in M101, may arise in HNRs. The upcoming 
AXAF and XMM will be able to spatially resolve such remnants and to obtain high 
spectral resolution data over the broad 0.1-10 keV range. The X-ray data, combined with 
detailed observations of the remnants in other wavelength bands, will provide important 
constraints on the rate, blastwave energy distribution, and progenitor type of the events. 
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Most excitingly, such a study may open up a new avenue to explore the origin and physics 
of GRBs. 

The author is indebted to Y.-H. Chu for her sharing insights and valuable information 
about the remnants during the course of this work, which is supported partly by NASA 
LTSA grant NAG5-6413. He is also grateful to A. Loeb, S. Oey, and B. Paczyhski for 
stimulating discussions on GRBs as well as on hypernovae and their remnants. 
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Table 1: X-ray-emitting Remnants Candidates in M101 



SNR 


R.A. 


Dec. 


Rg 


Ax — O 


r err (la) 


L(Ha) 


[SII] 


Size a 




Optical 


RHRI 


L x 


Name 


(h m s) 


(° ' ") 


(kpc) 


(") 


(") 


(ergs s -1 ) 


/Ha 


(PC) 


(km s -1 ) 


Shape 


Source 


(ergs s _1 ) 


MF37 


14 3 04.1 


54 27 36 


14 


3.3 


3.0 


4.3E37 


0.72 








H19 


6.9E38 


MF54 


14 3 20.7 


54 19 42 


3.6 


6.2 


3.1 


3.5E37 


0.52 


27 




filled 


H29 


1.4E38 


MF57 


14 3 24.2 


54 19 44 


4.2 


4.0 


3.1 


1.1E37 


0.47 


40 




filled 


H30 


1.1E38 


MF83 


14 3 35.9 


54 19 24 


7.7 


1.8 


3.1 


1.7E38 


0.76 


267 


~ 86 


shell 


H36 


1.2E38 


NGC5471B 


14 4 29.1 


54 23 53 


24 


3.1 


3.4 


4.0E38 6 


0.51 6 


60 b 


> 175 c 


shell 


H49 


3.0E38 



a The sizes of MF SNRs are recalculated, using the M101 distance of 7.2 Mpc instead of 

5.4 Mpc assumed by MF97. 

b from Chen & Chu (1999). 

c from Chu & Kennicutt (1986). 
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Fig. 1. — RHRI intensity contour map overlaid on an optical image of M101 from Palomar 
Sky Survey. The map is adaptively smoothed with a Gaussian of size adjusted to achieve 
a constant signal-to-noise ratio of ~ 6. The contours are at 1.1, 1.6, 2.4, 5.3, and 
12 x 10~ 3 counts s _1 arcmin" 2 ; the lowest level is about 3a above the local X-ray background. 
The locations of the X-ray-emitting remnant candidates (Table 1) are marked. 



Fig. 2. — 68%, 90%, and 99% confidence contours of the spectral parameters: kT vs. 
normalization of the Raymond & Smith thermal plasma fit to the X-ray spectral data of 
NGC5471B. The cross marks the best-fit parameter values. 
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